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We propose that future observation of supernova relic neutrino (SRN) can be used to probe 
neutrino decay models. We focus on invisible (e.g., Majoron) decays, and work out the general 
solution of SRN kinetic equations in the presence of oscillations plus decay. We then apply 
the general solution to specific decay scenario, and show that the predicted SRN event rate 
can span the whole range below the current experimental bound. Therefore, future SRN 
observations will surely have an impact on the neutrino decay parameter space. 

1 Introduction 

In general, massive neutrinos can not only mix, but also decay. The most stringent and safe limit 
comes from the nonobservation of decay effects on solar v flux. However, due to the relatively 
small distance from the Sun, this limit is very weak, namely Ti/rrn > 5 x 10~ 4 s/eV, where T{ 
represents the lifetime of the eigenstates Vi with mass rrii. Therefore, the possibility of neutrino 
decay with longer lifetimes (and coming from other astrophysical sources) cannot be excluded. 

Here we focus on the diffuse v background produced by all past core-collapse Supernovae 
(SN) in the Universe — the so-called supernova relic neutrinos (SRN). Future observations of 
SRN can probe decay lifetimes of cosmological interest. In fact, for SRN decay effects to be 
observable in our universe, it must be roughly TiE/rrii < 1/Hq, where Hq is the Hubble constant. 
By setting Hq = 70 kms -1 Mpc -1 and taking E ~ O(10) MeV (i.e., in the energy range probed 
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by supernova neutrinos), a rough upper bound for the "SRN neutrino decay observability" is 
obtained, namely Tj/m, < O(10 n ) s/eV. The comparison of this bound with the previous 
solar limit implies that SRN leave many decades in Ti/rrii open to experimental and theoretical 
investigations. 

Motivated by this challenging opportunity, in this talk we aim at showing how to incorporate 
the effects of both flavor transitions and decays in observable SRN spectra, by solving the SRN 
kinetic equations for two-body non-radiative v decay. This work is based on the results obtained 
in^, to which we refer the interested reader for further details. 

The plan of the talk is as follows. In Sec. 2 we discuss the general case of 3^ flavor transitions 
followed by decays, and give the explicit solution of the neutrino kinetic equations for generic 
decay parameters. Specific numerical examples (inspired by neutrino-Majoron decay models) 
are given in Sec. 3, in order to show representative SRN event rates and energy spectra in the 
presence of decay. Finally, in Sec. 4 we draw the conclusion of our work. 

2 Three-neutrino flavor transitions and decays 

In this Section we discuss and solve the neutrino kinetic equations in the general case of 3v 
flavor transitions plus decay. Notice that, for Ti/rrii values above the solar v bound, SRN flavor 
transitions occur in matter (and become incoherent) well before neutrino decay losses become 
significant, so that hypothetical interference effects between the two phenomena can be neglected. 
For our purposes, flavor transitions inside the supernova can thus be taken as decoupled from 
the subsequent (incoherent) propagation and decay of mass eigenstates in vacuum. 

2.1 3v flavor transitions 

We assume the active 2>v oscillation scenario. In this framework, the 2>v squared mass spectrum 
can be cast in the form (mf , mf, m|) = M 2 + (— ^j", +^f~, ±Am 2 ^ , where M 2 fixes the absolute 
v mass scale; Sm 2 and Am 2 govern two independent v oscillation frequencies, with |Am. 2 | 3> 5m 2 , 
as indicated by current data. The case of +Am 2 (—Am 2 ) characterizes the so-called normal 
(inverted) hierarchy. The elements of the mixing matrix U are parametrized in terms of three 
mixing angles (6>i 2 , #13, #23)- 

The yield a Y Vi of the i-th mass eigenstate at the surface of the supernova can be calculated 
by taking into account the stellar matter effects on v propagation in SN. These effects are 
parametrized in term of a level crossing probability Pjj = Ph (Am 2 , #13) among the instantaneous 
eigenstates of the Hamiltonian (the so called "matter eigenstates") in the dense medium b . The 
final results for the yields Y Vi at the exit from the supernova are collected in Table I in^. 

2.2 3v decays 

At the exit of SN, u$ mass eigenstates evolve indipendently until they reach the surface of the 
Earth. However, in their propagation in vacuum, they may decay. The number density of mass 
eigenstates n Vi (E, z) per unit of comoving volume and of energy at redshift z can be obtained 
through a direct integration of the neutrino kinetic equations, as described below. 
For ultrarelativistic relic neutrinos V{ the kinetic equations take the form 

n Ui {E,t)=R s ^{t)Y Vi {E)+ qji(E,t)-Ti^n Vi (E,t) . (1) 

m,j>m,i 

"The yield of a v species represents the time-integrated v luminosity. 

b In principle, there could be an additional level crossing probability Pl = PL(8m 2 , 612). However, according to 
current v phenomenology, it is Pl — for typical supernovae density profiles, and so we will neglect it hereafter. 



The left-hand side (l.h.s.) of Eq. (1) represents the Liouville operator for ultra-relativistic v, 
where H{t) is the Hubble constant at the time t. The right-hand side (r.h.s.) of Eq. (1) contains 
two source terms and one sink term. The first source term quantifies the standard (decay- 
independent) emission of V\ from core-collapse supernovae, which depends on the supernova 
formation rate Rsjsr(t)- The second source term, in which 
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dE'n Uj (E', t) B{ Vj -> Vi )Tj -± ^^(E', E) , (2) 

quantifies the population increase of U{ due to decays from heavier states Uj, with decay width 
Tj = 1/tj, branching ratio B(i/j — > i/j) and normalized decay energy spectrum i\) v .^ rVi . The last 
(sink) term on the r.h.s. of Eq. (1) represents the loss of Vi due to decay to lighter states with 
width Tj. 

Our result is that these equations can be directly integrated, by rewriting them in terms of 
the redshift variable z = z(t) and of a rescaled energy parameter e = e(E, z). By replacing back 
the variable E = e(l + z), one obtains the general solution of the neutrino kinetic equations, 
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where we have introduced the auxiliary function 
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In practice, these equations can be integrated numerically by following the decay sequence, 
i.e., starting from the heaviest state (qji = 0) and ending at the lightest state (Tj = 0). 

Finally, the v e flux at Earth (redshift z = 0), relevant for the inverse (5 decay reaction 
z/ e + p — > n + e + in a Cherenkov detector, is given by 

n De (E) ~ cos 2 On n„ x (E, 0) + sin 2 9 12 n P2 (E, 0) . (5) 



3 Applications to scenario inspired by Majoron models 

In this section we apply the general results of Sec. 2.2 to some representative decay scenarios, 
inspired by Majoron models. We consider only nonradiative (invisible) decays of the kind — ► 

{ Vj + X, where an heavier neutrino Vi decays into a ligher detectable (anti)neutrino plus an 
invisible massless (pseudo) scalar particle X, i.e. a "Majoron". 

We examine and compare a few representative 3v decay cases, which provide SRN yields 
higher, comparable, or lower than for no decay. For simplicity, we focus on two phenomeno- 
logically interesting cases in which the branching ratios and decay spectra become model- 
independent, namely, the case of quasidegenerate (QD) neutrino masses (m^ ~ nu > - rrij) 
and of strongly hierarchical (SH) neutrino masses (m, 3> rrij ~ 0). Table II in^ displays the 
relevant characteristics of the QD and SH cases. 



3.1 Three- family decays for normal hierarchy and quasidegenerate masses 

This decay scenario provides SRN densities generally higher than for no decay. The relevant fea- 
tures of this scenario are graphically shown in the left panel of Figure 1. The QD approximation 
forbids decays of neutrinos into antineutrinos and vice versa (see Table II iiM*). 
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Figure 1: Supernova relic v e spectrum (middle panel), and associated positron spectrum from v e + p — > n + e + 
(right panel), for a decay scenario with normal hierarchy and quasidegenerate masses (left panel, with r/m and 

branching ratios assignments). 

By construction, the decay scenario considered in this section is thus governed by just one 
free parameter (r/m). Notice that, for r/m ~ O(10 10 ) s/eV, SRN decay effects are expected 
to occur on a truly cosmological scale. For much larger values of r/m, the no-decay case is 
recovered. For much smaller values of r/m, SRN decay is instead complete, all SRN being in 
the lightest mass eigenstate v\ at the time of detection. 

Figure 1 shows the supernova relic v e energy spectrum, and the associated (observable) 
positron spectrum, for the considered decay scenario. The energy spectra for complete decay 
(red solid curves) appear to be a factor of ~ 2 higher than for no decay (black solid curve). 

For incomplete neutrino decay (i.e., for r/m ~ O(10 10 ) s/eV), one expects an intermediate 
situation leading to a SRN flux moderately higher than for no decay. Figure 1 displays the 
results for a representative case (r/m = 7 x 10 10 s/eV, red dotted curves). In conclusion, the 
decay scenario examined in this section can lead to an increase of the SRN rate, as compared 
with the case of no decay. The enhancement can be as large as a factor ~ 2, the larger the more 
complete is the decay. 



3.2 Three-family decays for normal hierarchy and m\ ~ 

This decay scenario provides observable SRN densities generally comparable to the no-decay 
case. In this case, the approximation of strong hierarchy (SH) can be applied to the decays of 
z^2,3 (and of ^2,3)- Figure 2 shows the supernova relic v e and positron spectra for this scenario 
where, as depicted in the left panel, all decay channels are open. This complex decay chain 
produces a substantial enhancement of the SRN energy spectrum at low energy, visible as a 
"pile-up" of decayed neutrinos with degraded energy in the middle panel of Fig. 2. Analogously, 
the case of incomplete decay (e.g., r/m = 7 x 10 10 s/eV, green dotted curves), is appreciably 
different from the cases of no decay and of complete decay only at low energy. 

In this scenario, the interesting effects of decay are almost completely confined to low v e 
energies, and are thus washed out in the observable e + spectrum, due to the cross section en- 
hancement of high-energy features. In conclusion, the e + spectra for the three cases of complete, 
incomplete, and no decay, turn out to be very similar to each other (right panel of Fig. 2), and 
so in this scanario decay effects cannot easily be distinguished through future SRN observations. 
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Figure 2: Supernova relic u e spectrum (middle panel), and associated positron spectrum from u e + p — > n + e + 
(right panel), for a decay scenario with normal hierarchy and mi ~ (left panel, with r/m and branching ratios 

assignments) . 



3. 3 Three-family decay for inverted hierarchy 

We conclude our survey of 2>u decays by discussing a scenario where the SRN density is generally 
suppressed, as compared with the case of no decay. The main features of this scenario are 
graphically shown in the left panel of Figure 3. It results that only the decay vi — > v\ (where 
the QD approximation is applicable) is relevant to SRN observations. In fact, decays to v% 
provide a negligible amount of z/ e 's (rx sin 2 613 < 1CP 2 ), so that the absolute value of mz makes 
no difference. In the above scenario, the case of complete decay (r/m <C O(10 10 ) s/eV) is trivial: 
since the final state is populated only by (and 1/3), the relic density of v e is negligibly small 
(rx sin 2 #13). The nontrivial case of incomplete decay (r/m ~ O(10 10 ) s/eV) is then expected 
to lead to an intermediate suppression of the SRN density. Figure 3 shows the numerical 
results for the specific value r/m = 7 x 10 10 s/eV, in both cases Pjj = 1 (solid curves) and 
Ph = (dashed curves). The neutrino spectra for incomplete decay (blue curves) appear to 
be systematically lower than the corresponding no-decay spectra (black curves), although the 
difference is mitigated in the positron spectra (right panel of Fig. 3). 

In conclusion, in the inverted hierarchy scenario of Fig. 3 the SRN signal is generally sup- 
pressed by neutrino decay, and can eventually disappear for complete decay. 



3.4 Overview and summary of 3v decay 

We think it useful to show also the behavior of the SRN signal for continuous values of the 
free parameter r/m, for the three different decay scenario discussed above. Figure 4 shows the 
positron event rate integrated in the energy window E e + G [10 , 20] MeV which might become 
accessible to future, low-background SRN searches^. For each scenario, the rate is normalized to 
the standard expectations for no decay and normal hierarchy (NH), and is plotted as a function of 
r/m. It results that neutrino decay can enlarge the reference no-decay predictions for observable 
positron rates by any factor / in the range ~ [0 , 2.3], depending on the particular decay scenario 
and provided that r/m is in the cosmologically interesting range below O(10 n ) s/eV. 

Since the current experimental upper bound on the SRN flux from SK is just a factor of 
~ 2-3 above typical no-decay expectations^, future observations below such bound are likely 
to have an impact on neutrino decay models. If experimental and theoretical uncertainties can 
be kept smaller than a factor of two (a nontrivial task), one should eventually be able to rule 
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Figure 3: Supernova relic v e spectrum (middle panel), and associated positron spectrum from v E +p — > n + e + 
(right panel), for a decay scenario with inverted hierarchy and generic 771,3 (left panel, with r/m and branching 

ratios assignments). 
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Figure 4: Positron event rates in the energy range [10 , 20] MeV for various decay scenarios, normalized to 
standard expectations for normal hierarchy and no decay. Notice how the expectations branch out (and then 
reach the complete decay limit) in the cosmologically relevant range r/m < 10 s/eV. 



out, at least, either the lowermost or the uppermost values in the range / £ [0 , 2.3], i.e., one of 
the extreme cases of "complete decay." Optimistically, one might then try to constrain specific 
decay models and lifetime-to-mass ratios through observations. 

4 Conclusion 

Neutrino decays with cosmologically relevant neutrino lifetimes [T;/mj < O(10 n ) s/eV] can, 
in principle, be probed through observations of supernova relic v e (SRN). We have shown how 
to incorporate the effects of both flavor transitions and decays in the calculation of the SRN 



density, by finding the general solution of the neutrino kinetic equations for generic two-body 
nonradiative decays. We have then applied such solution to three representative decay scenarios 
which lead to an observable SRN density larger, comparable, or smaller than for no decay. In 
the presence of decay, the expected range of the SRN rate is significantly enlarged (from zero 
up to the current upper bound). Future SRN observations can thus be expected to constrain 
at least some extreme decay scenarios and, in general, to test the likelihood of specific decay 
models, as compared with the no-decay case. 
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